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We develop an approach for self-consistent ac quantum transport in the presence of time- 
dependent potentials at non-transport terminals. We apply the approach to calculate the high- 
frequency characteristics of a nanotube transistor with the ac signal applied at the gate terminal. 
We show that the self-consistent feedback between the ac charge and potential is essential to prop- 
erly capture the transport properties of the system. In the on-state, this feedback leads to the 
excitation of plasmons, which appear as pronounced divergent peaks in the dynamic conductance at 
terahertz frequencies. In the off-state, these collective features vanish, and the conductance exhibits 
smooth oscillations, a signature of single-particle excitations. The proposed approach is general and 
will allow the study of the high-frequency characteristics of many other low-dimensional nanoscale 
materials such as nanowires and graphene-based systems, which are attractive for terahertz devices, 
including those that exploit plasmonic excitations. 
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I. INTRODUCTION 

A fundamental understanding of the physical processes 
controlling the complex space- and time-dependent be- 
havior of the carrier dynamics in reduced dimensions is 
essential to assess the technological potential of a variety 
of nanomaterials for future high-speed electronic devices. 
Such assessment, however, requires both experimental 
and theoretical techniques by which high-frequency ma- 
terial properties, such as the dynamic (ac) conductance, 
can be measured or calculated. 

Experimentally, much progress has been made in devel- 
oping techniques to probe the rf response of nanomateri- 
als. Among these nanomaterials, carbon nanotubes have 
received much attention due to their exceptional elec- 
tronic transport properties at dc^— and the hope that 
these carry over to high frequencies. Measurements of the 
high-frequency characteristics of carbon nanotube field- 
effect transistors (NTFETs)£~— indicate little decrease in 
performance up to gigahertz (GHz) frequencies. In addi- 
tion, the carrier dynamics in the terahertz (THz) regime 
was recently probed using time-domain techniques, sug- 
gesting that the carrier dynamics is determined by single- 
particle rather than plasmonic properties^ 

Theoretically, the problem of time-dependent trans- 
port has been approached using a variety of tech- 
niques such as scattering matrix theory^ - — Floquet 
methods)^— Boltzmann transport theory^ - — and 
nonequilibrium Green Functions (NEGF)^r— Even 
though the NEGF technique has become to some extent 
the standard in modeling electronic quantum transport, 
its application to time-dependent problems has been 
mainly focused on simplified few-level models ^ ~ 48 ' 51 " — 
or to few-atom one-dimensional wires or molecules i£Ir— 
While such model systems are invaluable to gain insight 
into the basic dynamic processes of simple quantum sys- 
tems, it can be difficult to relate them to more realistic 



devices for three main reasons. First, the self-consistency 
between the charge and the potential is needed to prop- 
erly determine the operation of the device in the presence 
of dynamic potentials. Second, it is necessary to cap- 
ture the rather complex dielectric environment of real 
devices. Third, most approaches have focused on ap- 
plying time-dependent signals at the source-drain, i.e. 
transport terminals ; 29 ' 36 " — rather than at the gate, a 
non-transport terminal. Physically, non-transport termi- 
nals do not supply the device region with charge, unlike 
source-drain contacts, but are coupled to the device chan- 
nel only through the (self-consistent) dynamic potential, 
which plays a similar role as the pumping potential in 
the field of parametric pumping] 39 : 49 

In this work, we make a first step towards solving 
this problem and develop a linear response theory for 
ac quantum transport employing nonequilibrium Green 
functions solved self-consistently with Poisson's equation, 
when a time-dependent signal is applied at the gate ter- 
minal. We apply the approach to a NTFET and de- 
termine its high-frequency response, showing that the 
approach cannot only describe time-dependent, single- 
particle quantum transport effects, but also is able to 
capture the plasmonic excitations of the device. 

The program of the paper is as follows: in Sec. II, 
we detail the formal theory for ac quantum transport 
and derive an effective Dyson equation describing the 
dynamics of the system for a time-harmonic signal at a 
non-transport terminal. Special attention is given to the 
practical calculation of the frequency-dependent charge 
density for which we develop a computationally efficient 
scheme, a prerequisite for calculating the self-consistent 
ac response of larger systems, as we have demonstrated 
previously^ In Sec. Ill, general expressions for the ac 
particle current and associated conductance are derived 
along with a brief outline of the current partitioning 
scheme^ and how it applies to a multi-terminal device 
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with non-transport terminals. In Sec. IV, we apply the 
theory to a NTFET. There, we discuss details of the sig- 
nificance of the operation mode of the device, and the 
self-consistent feedback between charge and potential for 
collective excitations. Our conclusions are presented in 
Sec. V. 



II. GENERAL APPROACH 

In this section we describe the development of the ac 
approach, which consists of three steps: 1) definition of 
the model Hamiltonian of the total system, 2) formulat- 
ing the quantum dynamics and non-equilibrium statistics 
in terms of Green functions in the energy domain, and 3) 
self-consistent calculation of the ac charge and potential. 

A. Model Hamiltonian 

We begin by specifying the Hamiltonian operator of 
the system. As usual, the total system is divided into 
three isolated regions following the partitioning scheme of 
Caroli and co-worker a 56 ' 57 The Hamiltonian of the entire 
infinite system is written as 

H = H d + H c + H t , (1) 

where Hd is the Hamiltonian for the device region, H c 
refers to the two semi-infinite leads, and H t couples the 
device region to the leads. In a site representation the 
device Hamiltonian is given by 

H d = H°+H? c + Ht c , (2) 

where 

H d = 6° & n C n + — 2^ ti^mC^Cm + h.c. , (3) 

n n,rn 

and 

H? C = J2 U n C &n , Hf C =J2Un(t)cic n , (4) 

n n 

where c\ and c n refer to fermionic creation and annihi- 
lation operators at site n. defines the equilibrium 
electronic structure of the isolated system. The electron- 
electron interaction is approximated on the Hartree level 
and has two components H® and H^ c ', cf. Eq ((2J). 
The term U® c represents a spatially-varying, but time- 
independent electrostatic potential, such as the one 
present when calculating the DC properties and leads 
to a rcnormalization of the onsite energies e° . 

The new physics studied here originates from the pres- 
ence of an a priori unknown time- and space- dependent 
potential U n (t) induced by externally applied time- 
dependent fields. As further discussed below, both U® c 
and U n (t) must be determined separately by solving Pois- 
son's equation in a self-consistent manner. In general, the 



approach allows to investigate the dynamic response be- 
yond the Hartree approximation of the Coulomb interac- 
tion by including exchange and correlation functional^ 
calculated self-consistently. 

The Hamiltonian for the two contacts to the left and 
right (a = s, d) of the device reads 

H c = Yl e k a c\ a c ka , (5) 

k , OL 

where c ka and Ck a are fermionic creation and anihilation 
operators for a particle in terminal a in state k. We note 
that equations ([3]), (@|, and ([3]) differ from those con- 
sidered previously where a time-dependent source-drain 
bias is considered, in which case the onsite energy e° Q 
of the contacts become time-dependent rather than the 
ones of the device. 

Finally, the Hamiltonian 

H t=Yl T n,k a c\c ka + T* M c{ a c n , (6) 

ka,n 

couples the device sub-space with the semi-infinite source 
and drain reservoirs, and allows for a physical exchange 
of particles through the device-contact interface. There- 
fore, the tunneling Hamiltonian Eq. ((BJ) describes only 
the coupling between the device and transport terminals, 
but not to non-transport terminals. 



B. Quantum dynamics and nonequilibrium 
statistics 

The next step is to describe the carrier dynamics 
within the device scattering region using Green func- 
tions. The Green functions arc in general functions of 
both space and time, e.g. G(rt;r't'). However, to sim- 
plify the equations for compactness we adopt a short- 
hand notation G(t,t') = G(rt;r't'). In addition, when- 
ever regular functions appear with Green functions in the 
same equation, we also omit the spatial dependence on 
the regular functions. 

We start with the time-dependent Dyson equation^ 

C{t,t') = gZ{t,H) (7) 
+ Jdt 1 dt 2 9l{tM)^ '(ti,t 2 )G~< '(ta.O , 

where g^(t,t') = g^{t — t') refers to the 
retarded/advanced (7 = r,a) Green function of 
the isolated system. The self-energy E 7 accounts 
for all interactions of the isolated system with its 
environment. In our case, the self-energy £ 7 (£, t') can 
be divided into three contributions 

£ 7 (M')= K(t-t') + U DC S(t-t') + U(t)6(t-t') . 

a — s,d 

(8) 
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The first term S c = J2 a =s d ^ a * s the contact self- 
energy and corresponds to the quantum-transport open- 
boundary conditions connecting the device region with 
the semi-infinite source and drain contacts. The sec- 
ond term is a scalar potential and represents the inter- 
nal response of the device to externally applied time- 
independent fields. The third term is the prominent 
feature in the ac theory presented here, and describes 
the dynamic response of the device due to external time- 
dependent fields. Contrary to most studies where the ac 
signal is applied at the source-drain terminals ; 32 ' 47 ' 50 ' 54 
in our case the time-dependent signal is applied at the 
gate terminal. This implies that the induced potential 
U(t) distorts only the device scattering region, while the 
contacts remain in steady-state. 

We now switch from the time-domain into energy- 
representation through a double-time Fourier-transform 
defined as^ 7 - 

FiE.E^f^'e^e-'^FM (9) 

and 

F(t,t')= f ^-^- e - iEt l n e iE,t '> n F{E,E') , (10) 

so that the self-energy, cf. Eq. is given by 

^(E, E') = 2ttS(E - E') [T%(E) + U DC ] + U(E - E'). 

It is worthwhile mentioning that in energy domain the 
contact self-energies are local in energy, reflecting that 
under steady-state conditions there is no mixing between 
states with different energy within the reservoirs. On the 
other hand, the original time-local potential U becomes 
now in energy domain non-local, implying that a time- 
dependent potential mediates transitions between states 
at different energies within the device scattering region. 

Fourier transforming Eq. (O and using (|TTj) . one de- 
rives an effective Dyson equation for the device 

G^^E.E 1 ) = 2ttS(E - E')Gl{E) (12) 

+ / ^Gl{E)U{E-E)G^E,E') , 

where 

Gl{E) = [g^E)- 1 -U DC -^(E)]- 1 , (13) 

and 

92(E) = [(E±i V )I-H$\- 1 , (14) 

with an infinitesimal rj > 0. What we have gained in 
re-formulating Dyson's equation is to partition the full 
dynamic response of the system described through the 
two-energy Green function G(E, E') into its dc and ac 
components given by the first and second term in equa- 
tion (|12[) . respectively. Importantly, the DC component 



determined by the newly defined Green function Gq, cf. 
Eq. (fT3"|) . refers no longer to the response of the isolated 
system <?q , cf. Eq. (|14p . but rather describes the system's 
response in contact with the leads and subject to a dc 
electrostatic potential. Hence, Gq defines the operation 
point of the open system under DC steady-state. The ac 
component, i.e. the second term in Eq. (|12j) contains this 
term as well and determines the distortion of the system 
away from the operation point Gq, and is driven by the 
time-dependent potential U(t) leading to a coupling of 
states at different energies. 

We still need to know how the total nonequilibrium 
particle distribution G < deviates from its (reference) dis- 
tribution at dc in the presence of the ac potential U . 
This is accomplished by mapping Dyson's equation for 
G < , symbolically written as G < = [Go + GoUG] < , onto 
the real-time axis utilizing the Langreth rule s 58 ' 59 which 
gives: G< = G< +G<UG a + G r UG< . This integral equa- 
tion can be solved exactly making use of Eqs. (|T2"j) and 
([TB")) . Details of the derivation are found in appendix A. 
After Fourier transform the particle distribution is given 
by 

G<{E,E') = 2irG${E)5(E - E') (15) 
+ J^[G<(E)U(E-E)G a (E,E') 

+ G r (E,E)U(E~E')G<(E')] 

f dE x dE 2 dE 3 

+ / x 

J 2tt 2tt 2tt 

G r (E,E 1 )U(E 1 -E 2 )G<(E 2 ) x 

U(E 2 -E 3 )G a (E 3 ,E') , 

where Gf(E) = G r Q (E)T,f(E)Gf ) (E) corresponds to the 
nonequilibrium spectral particle density at dc. The func- 
tion £<(£) = J2 a = s , d iU(E)r a (E) where T a (E) = 
i(Y7 a — S°) is the broadening function, and f a (E) 
i 1 ■ e (£-Mo)/fci3Tj j g ^ e p erm j f U nction at tempera- 
ture T with \x a being the chemical potential of terminal 
a. 

While the set of equations (|T2"j) - (fl"5"|) developed so far 
describe entirely the quantum transport and nonequilib- 
rium statistics, they do not allow to determine the dy- 
namic potential U. This must be obtained by solving 
Poisson's equation 

V [e(r)VC/(r, E - E')} = -p(r, E - E') , (16) 

with the frequency-dependent charge density 

/j 771 
— G<(E+,E). (17) 

The calculation of the ac charge density using Eq. (fT5|) 
requires G < to be evaluated at two energies (E + ,E) = 
(E + fiu,E), in contrast to the dc case where only one 
energy is needed. 

Equations ()16[) and (jTTJ) implement the self- consistent 
coupling between electrostatics and transport, which rep- 
resents the key component in our ac approach. Note 
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that at the frequencies considered here the electromag- 
netic fields respond instantaneously, so that the full time 
dependence in Maxwell's equations can be neglected. 
Poisson's equation is supplemented by boundary condi- 
tions appropriate for the problem at hand, and e(r) is 
a space-dependent dielectric constant that can account 
for more complex inhomogcneous dielectric environments 
quite common in devices. 

The set of Eqs. ([12")) - ([TT)) describe the noncquilibrium 
quantum dynamics and its coupling to Poisson's equa- 
tion for an arbitrary time-dependent potential U, and 
can thus describe situations beyond linear response, in 
general. However, the numerical implementation of the 
full non-linear theory requires the calculation of a triple 
energy integral in Eq. (|15p . which is prohibitive at this 
time given the need for self-consistency to capture the 
plasmonic response of real devices as discussed in Sec. 
IV. 



C. Linearized equations 



To proceed further, we now apply a time-harmonic sig- 
nal at the gate terminal v g (t) = vq cos(wt) of small ampli- 
tude vo and frequency ui, and seek the potential response 
in the form U(r,t) = V(r,u>) cos(wi), which reads in en- 
ergy domain 

U(E) = i V(r, u>) [5{E + huo) + S(E - hu))] . (18) 

Keeping only terms to linear order in V, the ac transport- 
Poisson equations take the form 

C(E+,E) = 2vGl{E)5{hu) + \cP {E + )V{w)G%{E), 

(19) 

G<(E+,E) = 2TrG<(E)5(huj) + ^G<(E+)V(uj)G a (E) 



-G r (E+)V(u)G<(E) , 



(20) 



PM = le /^ G<(E+,E), (21) 
-p(r,u) = V[e(r)W(r,w)] . (22) 

D. Numerical calculation of p(r, w) 

An integral part in the self-consistent transport calcu- 
lations is the determination of the charge density. In 
practice, one has to evaluate the integral in Eq. (|2ip 
which is often performed by direct integration along the 
real energy axis. In many cases, this is a sufficient ap- 
proach because the spectral density of states has a fi- 
nite bandwidth, thus narrowing the integration window. 
However, such conditions are rarely realized in more real- 
istic device models. For instance, even in a simple tight- 
binding representation of a NTFET (see Sec. IV) the 



bandwidth of the valence and conduction band is about 
10 eV, in which case the calculation of the charge density 
through a real-axis integration can become prohibitive 
for self-consistent calculations even at dc. This becomes 
an even more severe bottleneck in the case of ac simu- 
lations, where now the charge has to be determined at 
every frequency uj. 

In the following, we describe a computational efficient 
approach, which permits the calculation of the frequency- 
dependent charge density p(ui) by exploiting contour inte- 
gration in the complex energy planei 60 ' 61 The basic idea 
is similar to the dc case, i.e. to separate in Eq. (|2~lj) 
the zero-bias contribution to p(ui) from the non-zero-bias 
component. If we further assume that the lowest chem- 
ical potential is at the drain terminal, i.e. p^ < pi 
the frequency-dependent particle distribution at zero- 
bias (ZB) reads 



G< B (E+,E) = -A+V(w)G%f+ + -G r >+V(os)A f d 



(23) 

where a + superscript indicates a function evaluated at 
E + %U3, and the absence of such a superscript indicates a 
function evaluated at E. The steady-state spectral den- 
sity is given by 



A 



(+) 



i G, 



G 



»,(+) 



(24) 



Gzb contains Fermi functions evaluated a two different 
energies E and E + , reflecting the nonequilibrium nature 
of the ac charge density for finite frequencies even at zero- 
bias, which means that an externally applied ac signal 
acts as if a frequency-dependent bias were applied. 

Taking advantage of this ac-signal-bias analogy and 
noting that f£ < fd, one can split again the zero- bias 
particle density into its equilibrium and noncquilibrium 
components, i.e. G^ B = G"^g q + G"^g eq , which after 
re-arrangement take the form 



Gfs q = - g {fd + fl) Im {G% + V(w)Gr } , (25) 



and 



Gfr = --{U- ft)[^{G^ + V{u)Gl\ (26) 

-G^V^G^} . 

The equilibrium part is analytic in the upper com- 

plex plane, since it consists of the product of two retarded 
Green functions Gq and Gq + each of which has poles 
only in the lower complex planed Therefore, the equi- 
librium zero-bias ac particle density, which involves all 
states below the frequency-dependent chemical potential 
Pa , can be efficiently calculated through integration over 
a complex energy contour . 60 ' 61 Conversely, the nonequi- 
librium component at zero-bias G^'g Lq is non-analytic, 
because both the retarded and the advanced Green func- 
tions are needed with their corresponding poles located 
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in the lower and upper complex plane, respectively. How- 
ever, this does not pose a serious problem in practice as 
the integration range is limited to a finite energy window 
given by Hoj, i.e. the difference between the chemical 
potentials na — A* d ~ 



III. AC RESPONSE FUNCTIONS: CURRENT 
AND CONDUCTANCE 

The set of equations (fTi?)) - ([2"2"j) developed in the pre- 
vious section allow the determination of the frequency- 
dependent Green functions, which can now be used to 
obtain ac response functions. One basic response func- 
tion to characterize transport is the dynamic conduc- 
tance <7ct/3, which relates the total ac current I a with the 
voltage applied at terminal /3. Under time-dependent 
conditions this conductance is not entirely determined 
by the particle current, but has in general contributions 
from the displacement current as well. In the following 
subsections, we derive an expression for the particle con- 
ductance, and summarize how displacement currents can 
be included in the total conductance. 



A. Particle current 

The first contribution to the total current consists of 
the flow of charged particles through the terminal a, and 
is hence determined by the dynamic change of the particle 
density at this terminal 

m = -4<#„(f)> = -e-§5>L • (27) 



Making use of the fermionic anti-commutator relations,— 
and the Heisenberg equation of motion for operators 



J 



O = j- [H, O] with H the total system Hamiltonian, cf. 
Eq. (TT|), one derives the well-known matrix equation for 
the equal-time particle curren t 27 ' 46 ' 50 

P a {t) = \Tr J dt> [G r (t,t')X<(t>,t) -X<(t,t')G a (t',t) 
+ G<(t,t')X a a (t>,t)-X r a (t,t>)G<(t',t)] . (28) 
Its corresponding energy representation read s 27 ' 46 ' 47 
dEdE' 



H Tr J ^[ G< ( E ^ E ')K(E',E) 

^ a (E+,E')G < (E\ E) + G r {E+,E')T,<{E',E) 
Z<(E + ,E')G a (E',E)] . (29) 



Equation (|29[) simplifies further if we exploit the steady- 
state property of the contact self-energies from Eq. (fTTj) in 
which case one obtains the frequency-dependent particle 
current 

P a (u) = ~Tr J dE[G<(E+,E)j: a a (E) 

- ^ r a {E + )G < {E + ,E)+G r {E+,E)^<{E) 
- Z<(E+)G a (E+,E)] . (30) 

We note that the expression for I^(uj) differs from 
those derived in Refs i 46 ' 47 ' 50 , since those applied a time- 
dependent voltage at the source-drain, instead of the gate 
excitation considered here. 

One can now derive the dynamic particle conductance 
by expanding and G < appearing in Eq. (|30|) to linear 
order in the terminal voltage vp, and utilizing Eq. (|19p 
and (2D) to substitute for G^{E+,E) and G < (E+,E). 
These linearized expressions are summarized in Appendix 
B. Inserting all relevant terms in Eq. ([30]) and keeping 
components linear in V(uj), we derive for the frequency- 
dependent particle current 



= YhH Tr [ dE^+V^G^ -^ + G a + V(u;)G a }5 aP 

+ G< : +V{u)G^ a + G r o ' + V(u,)G< E a a - Z%+G<£V(w)G5 - ^ a + G r , + V(u)G< p 

I 



(31) 



By definition, the (tensor) prefactor that relates the ter- response particle conductance, and can be read directly 
minal current I a with the applied bias vp is the ac linear from Eq. (|31[) : 



!±Tr 
2 h 



dE 



'{gT ' + V(uj)G^< - Z<' + G a + V(u)G a } S a p 
+ G<;+V(u)G^ a a + G r + V(Lo)G<^ a a - E r ^ + G<;+V (u;)G a - ^ a +G r Q + V^)G< fi 



(32) 



6 



B. Displacement current I d ((jj) 

Under time-dependent conditions the particle con- 
ductance does not in general obey sum-rules, i.e. 
J2 a 9aP = an d J2/3 9aP = 0, reflecting current con- 
tinuity and gauge-invariance, because the displacement 
current present under ac conditions is often discarded. 
The current partitioning scheme of Wang et al£l allows 
to re-establish these sum-rules by taking displacement 
currents into account. 

The basic idea of this scheme can be summarized as 
follows: starting from the charge continuity equation, 
d t p + V • j p = 0, and integrating over the volume one 
obtains Kirchoff's current law, i.e. I (t) + J2 a ^a = 0- 
I? refers to the particle current through terminal a, and 
can be associated with a particle conductance through 
1% = ^2pg a pVp with vp- the voltage at terminal /3. The 
displacement current I d (i) = dtQ{i) accounts for the dy- 
namic change of the total charge, and is non-zero under 
time-dependent conditions. 

To obtain an expression for the total conductance de- 
fined by I a = 9ap~vp- one needs to know how the cur- 
rent I a is split between the particle and the displacement 
current at each terminal. While the particle component 
IP is directly accessible through transport, this is not im- 
mediately possible for I d , since only the total rather than 
the terminal displacement current is known. This prob- 
lem can be resolved by making two Ansatze for the termi- 
nal and total displacement current,— i.e. I a := IP+A a I d 
and I d := ^ ^ g^v p , where g d defines the displacement 
conductance, and permits to specify a total conductance: 
9a/3 = 9 a fj + A a g d . The partitioning factor A a can be 
determined by employing the sum-rules ^2 a g a /3 = and 
Sfl 9<x& ~ 0; so th & t * ne total conductance is given by^£ 



9*0 = 9a.fi 



S 7 9aj j 

Zd9f} 



(33) 



7 a 7 



and constitutes a (TV x TV)-matrix for a system with TV- 
terminals, in general. 



IV. APPLICATION: NANOTUBE FET 

In this section, we apply the approach developed in 
the previous sections to a ballistic nanotube transistor 
shown in Fig. Q] with a channel length of L = 20 nm. 
In general, the high-frequency properties of this three- 
terminal device can be determined by any component 
of the (3 x 3)-conductance matrix, cf. Eq. (|3"3")l : here 
we chose the source-drain conductance g s d(u) calculated 
at zero-bias and T = 300 K. Since the ac signal is ap- 
plied only at the gate terminal, which couples capac- 
itively to the nanotube channel, one can assume that 
the dominant contribution to the displacement current 
is carried by the gate, i.e. I d « I g . In this case, the 
partitioning factors for the three terminals simplify, i.e. 




FIG. 1: (Color online). Cross section of the NTFET device 
in cylindrical geometry and embedded in a dielectric medium. 
The dashed rectangle specifies the computational domain. 



(A s , Ad, A g ) « (0,0, 1), so that the source-drain ac con- 
ductance is given by g s d(^>) ~ 9sd( w )- Here, we focus on 
one particular channel length to discuss aspects of the 
methodology that are essential for the proper descrip- 
tion of the ac behavior. The properties of such devices 
with different dimensions were presented by us in detail 
in Ref^. 



A. Transistor response in the dc operation point 

We begin our analyis by specifying the setup of the de- 
vice shown in FigurcQ]following the modeling approach of 
RefjSs. The channel, which consist of a semi-conducting 
tube with chirality (to, n) = (17, 0) and radius 0.66 nm, is 
placed in the center of a cylindrical hole with radius 0.96 
nm and surrounded by a dielectric with a permittivity of 
16 (Hf0 2 ). 

The equilibrium electronic structure of the nanotube 
is described within a p z tight-binding model with diag- 
onal matrix elements e° !; = 0, and off-diagonal elements 

*2i,2i-l = *2*-l,2i = 27 cos (^); *2t,2t+l = *2i+l,2t = 7j 

where to refers to the number of carbon atoms per ring. 
The periodic boundary conditions along the tube circum- 
ference leads to a quantization of the wavefunction, so 
that the NT electronic structure can be classified by an 
angular momentum J = 1, . . . , m labeling the subbands. 
We chose 7 = 2.5 eV for the it carbon-carbon bond en- 
ergy, so that the bandgap between the highest valence 
and lowest conduction band (J = 6) is E g = 0.55 eVi^ 

The contacts are semi-infinite extensions of the NT 
channel, and described through self-energies Y7 a = ~f 2 g r a 
for each contact (a = s, d) where 7 couples the first/last 
ring of the NT channel to the surface of the contacts to 
the left and right.— The surface Green function g r a is cal- 
culated numerically at each energy using a matrix itera- 
tive schemed The matrix elements of the retarded Green 
function Gq for the NT channel are obtained employing 
a recursive algorithm^ The function of the two embed- 
ding metallic regions is to electrostatically dope the ends 
of the NT channel. In all simulations, the equilibrium 
Fermi level of the semi-infinite NT source/drain contacts 
Ep is set at —1.0 eV below the NT midgap energy before 
self-consistency, which gives after self-consistency p-type 
Ohmic contacts. 

Due to the cylindrical symmetry, the 3D Poisson's 
equation Eq. (|22p reduces to a two-dimensional (2D) 
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problem. In this case, Poisson's equation is discretized 
along the axial and radial axis within the 2D simulation 
domain as marked by the rectangular box using finite- 
differences^ and the resulting linear matrix system is 
solved by successive overrelaxation^ Along the domain 
boundary we impose homogeneous von Neumann bound- 
ary conditions for the electrostatic potential (VV = 0), 
and use Dirichlet boundary conditions (V = const.) at 
the perfect-metal source, drain, and gate terminals. Pois- 
son's equation requires a 3D charge density in real-space 
as input. However, an orthogonal tight-binding represen- 
tation of the NEGF transport equations calculates the 
total charge per NT ring. A 3D charge density can be 
obtained by smearing of the total charge per ring along 
the axial and radial direction of the 2D domain using 
Gaussian smearing functions^ 

The first step in determining the AC response of the 
transistor is to choose an operation point either in the off 
or on state, which is controlled by an appropriate DC gate 
bias V g . Figure [5] (a) shows the output characteristics for 
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FIG. 2: (Color online), (a) dc output characteristics for L — 
20 nm NTFET. The arrows mark the operation point in the 
on and off state, while panel (b) shows the respective self- 
consistent bandprofiles. The dashed-dotted line is the Fermi 
level with Ef = — 1 eV. 

our NTFET specified by the dc source-drain conductance 
G s d in the absence of a gate perturbation (vo = 0) with 
arrows marking the selected on and off states. At zero 
frequency, the conductance g s d{uJ = 0) examined in the 
forthcoming sections is related to G s d by its slope taken 
at the operation point, i.e. g s d(uJ = 0) ~ *gy d ^o- Fig - 
urc [21(b) displays the respective dc bandprofiles with the 
band being flat in the on state leading to a maximum con- 
ductance of 2e 2 /h (per spin) whereas in the off state the 
hole current is reduced due to the gate-controlled barrier 
in the channel. 
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FIG. 3: (Color online), ac response for a L = 20 nm NTFET 
in the off state: (a) real/imaginary part of the frequency- 
dependent conductance, and (b) color plot of the NTFET 2D 
density-of-states illustrating the resonant photoexcitation of 
carriers between spatially and energetically oscillating states. 
The solid black line marks the valence band edge, and the 
Fermi level is Ef = — 1 eV. 



due to the quantum interference of states by the barri- 
ers. Photocxcitations of carriers between states associ- 
ated with maxima in the DOS lead to maxima in g s d{uj), 
while its minima are caused by transitions between max- 
ima and minimal An oscillatory behavior of the conduc- 
tance is hence a signature of single-particle excitations, 
and is preserved when the self-consistent feedback is dis- 
abled as will be shown further below. We note that at low 
frequencies the real part of the conductance is negative. 
This is because in the limit uj — > the ac signal perturba- 
tion v g (t) = v cos(wi) becomes effectively a positive DC 
bias v g (t) = v > superposed to V g . According to the 
dc transfer characteristics of Fig. [2] (a) dG s d/dV g < 0, so 
that an increase of V g by vq leads to a reduction in the 
conductance. 



B. Transistor response in the off-state 

We now superpose an ac signal of small amplitude 
vo = 10 mcV and frequency uj to the dc gate bias V g . Fig- 
ure 0(a) shows the dynamic conductance in the off state, 
with real and imaginary parts having oscillatory charac- 
ter as a function of frequency uj. One can understand 
this behavior from the space- and energy-dependent 2D 
density-of-states (DOS) shown in Fig|3](b). For a given 
position z along the tube the DOS oscillates in energy 



C. Transistor response in the on-state 

The dynamic response is quite different in the on state 
as shown in Figure 0] (a). For small frequencies the dy- 
namic conductance is slightly negative for the same rea- 
son as in the off state, and exhibits a pronounced diver- 
gence at a discrete frequency of about rj 36 THz. Away 
from this resonance the conductance is oscillatory simi- 
lar to the off state, as shown more clearly in the inset of 
Fig. 11(a). 

In order to identify the nature of this resonance, we de- 
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FIG. 4: (Color online), ac response for a NTFET with L = 20 
nm in the on state. Panel (a) shows the real/imaginary part 
of g s d(w) with an enlarged view of the data shown in the inset. 
Panel (b): Amplitude of the normalized potential | V(oj) / Viait | 
with the dielectric function e{ui) shown in the inset. The 
arrows mark the near-resonance frequencies ui^ . Panel (c) 
shows the large amplitude potential profile V(lu) for . 



termine the response of the electrostatic potential V(lu) 
for fine-sampled frequencies near the divergent behav- 
ior. In Figure H] (b) we show the ratio of V(ui)/V ext 
with V ext the external perturbing potential. Interest- 
ingly, upon approaching the resonance the amplitude of 
the potential diverges, cf. Fig. |U(b). Alternatively, one 
can evaluate the frequency-dependent dielectric screen- 
ing e(uj) = V ex t/V(u>) shown in the inset of Fig. @] (b), 
which has a clear zero crossing at wo, while the poten- 
tial undergoes a change in sign as displayed in Fig. @] 
(c). These observations verify that the divergent behav- 
ior of the dynamic conductance observed in the on state 
is attributed to the excitation of plasmons which co-exist 
with the single-particle excitations. 



D. Importance of self-consistency 

In the previous section, we were able to identify the 
basic features in the dynamic conductance such as the 
oscillatory and divergent characteristics with the single- 
particle and collective behavior of the channel electrons, 
and concluded that plasmons can only be excited if the 
device is operated in the on state. 
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FIG. 5: (Color online). Comparison of the self-consistent 
(SC) vs. Non-self-consistent (Non-SC) response for a L = 20 
nm NTFET. Panels (a)-(c) is the response in the off-state, 
with panel (a) the conductance and panels (b) and (c) the 
charge density. Panels (d)-(e) show the behavior in the on- 
state. 

Are there other, more fundamental prerequisites irre- 
spective of the operation point, which determine whether 
the system can be driven into a collective state at 
all? The answer is yes, and is related to the self- 
consistency between charge and potential. In Figure [5] 
we compare the conductance g s d{w) and the frequency- 
dependent charge density p(z, w) calculated using the full 
self-consistent (SC) and a non-self-consistent approach 
where in the latter case the dynamic conductance is cal- 
culated in one step from the dc band profile. 

In the off state, the most apparent difference between 
the SC vs non-SC case is that, while the SC amplitude of 
the conductance is larger, the smooth oscillatory behav- 
ior is preserved as shown in Fig. [5] (a). Hence, the single- 
particle excitation spectrum is - at least qualitatively - 
not affected by the charge-potential feedback. This is also 
apparent in the ac charge density p(z,uj), cf. Figs. [5] (b) 
and (c), which exhibits space- and frequency-dependent 
oscillations in both the SC- and non-SC case. 
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In the on state, eliminating the feedback loop has 
a quite different impact on the response as demon- 
strated in Fig. [5](d). The plasmonic component visible 
through a divergent conductance g s d(u) with SC van- 
ishes for the non-SC calculation. This drastic change in 
the response from the (divergent) plasmon-dominated to 
single-particle characteristics is again clearly reflected in 
the ac charge density p(z,u>) calculated with (SC) and 
without (Non-SC) feedback shown in Figs. [5] (e) and (f). 
In the SC case, the charge density has a large amplitude 
at resonance with a peak in the middle of the channel, a 
feature that is absent in the non-SC calculation. 



V. CONCLUSIONS 

We develop an approach for ac quantum trans- 
port within the nonequilibrium Green function formal- 
ism, which allows to determine the frequency-dependent 
charge and potential under excitation at a non-transport 
terminal within a fully self-consistent framework. 

The capability of our approach to determine the high- 
frequency properties of systems in complex environments 
is demonstrated using a nanotube transistor with an ac 
signal applied at its gate terminal. In the off state, the 
dynamic conductance shows oscillations that originate 
from single-particle excitations between quantized energy 
levels. When the device is operated in the on state, the 
dynamic conductance exhibits discrete divergent peaks 
at terahertz frequencies. These peaks are associated with 
plasmonic excitations of the charge density at the reso- 
nant frequencies of the transistor acting as a quantum 
cavity. It is shown that the self-consistent coupling be- 
tween charge and potential is an essential component in 
the ac transport theory to capture plasmon excitations 
of the system. A non-self-consistent approach misses 
this important physics, and can only provide information 
about the single-particle excitation spectrum. 

The proposed approach is not limited to study the ac 
response of nanotube devices, but can be applied to ex- 
plore nonequilibrium, time-dependent electronic and op- 
tical processes in other low-dimensional materials such 
as nanowires, graphene, or molecules, including the ex- 
ploration of their collective excitation modes for novel 
plasmon-based nanoscale devices. 
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Appendix A: DERIVATION OF THE 
NONEQUILIBRIUM PARTICLE DENSITY G< 

In the following we detail the derivation for the parti- 
cle density G < (E,E'), cf. Eq. ([15]). We start from the 
expression for the (time-domain) Dyson equation for G < 
mapped onto the real axis utilizing Langreth rules ; 58 ' 59 
and symbolically written as 

G< = (Go + G C/G) < = G^ + GqUG ci + G r UG < (Al) 

where we have used that U K = U(t)S < (t — t')=0 for a 
time-local potential^ This equation can be re-arranged 
by collecting the Gq terms first 

G< = Gjf (1 + UG a ) + G r UG < . (A2) 

Equation (IA2[) can be solved through iteration by insert- 
ing the expression for G < on the l.h.s. into the second 
term on the r.h.s., and collecting now the Gq (1 + UG a ) 
elements. After the first iteration one obtains 

G< = (1 + G r U) G£ (1 + UG a ) + G r UG r UG < . (A3) 

The retarded Green function Gq in the prefactor (1 + 
GqII) is the first term in Dyson's scries, G r = G 7 Q + 
GqIIG 7 ' , which becomes more transparent when iterating 
one more time 

G< = (1 + G 7 Q U + G r UG , U)G< (l + UG a ) (A4) 
+ G r UG r UG r UG < 

= [l + (Gl + GlUG^)U]G<{l + UG a ) (A5) 

+ G r Q UG r UG r UG< . (A6) 

Iterating to infinite order, this Dyson scries converges 
towards G r , so that the final expression for the non- 
equilibrium particle density reads 

G< = (I + G r U) G< (1 + UG a ) . (A7) 

Equation ([15")) in section II. B corresponds to Eq. (|A7[) 
after Fourier transform. 



Appendix B: Small Bias Expressions for G < and E^f 

The conductance associated with the particle and dis- 
placement current are response functions which relate the 
terminal current I a with the terminal voltage vp in a lin- 
ear manner. In order to derive a formula for the con- 
ductance g p a p given in section III, one needs linearized 
expressions for S < and G < . These can be easily ob- 
tained from the Taylor expansion of the Fermi function 
fp = fp(E) = 1/ [l + e (E -^- 0+ev ^/ kBT ] to first order 
in the terminal voltages vp, i.e. 

fp ~ fp,o + efp vp , fj3 = - jr^ f/3,0 (1 - fp,o) , (Bl) 
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where np^ ls the chemical potential of terminal (3 at zero 
bias, and fp fi = 1/ [l + e( £ ^wl/ tsT '] the correspond- 
ing Fermi function. Inserting Eq. (|B1[) into and G < 
one obtains the following set of linearized expressions 

S| = ifpTp = £< + £< ^ (B2) 
s |,o = ir >//3,o , S< = il^/p (B3) 

and 



— ^O^^O^O i u 0,p — ^0^/3 U ' l co i 

j3=s,d 



with Tp = i[ So — So J the broadening function. 
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